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Rough surface effect on tlie Meissner diamagnetic current in tlie normal layer of proximity 
contact N-S bi-layer is investigated in the clean limit. The diamagnetic current and the screening 
length are calculated by use of quasi-classical Green's function. We show that the surface 
roughness has a sizable effect, even when a normal layer width is large compared with the 
coherence length ^ = vf/ttTc- The effect is as large as that of the impurity scattering and also 
as that of the finite reflection at the N-S interface. 
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§1. Introduction 

A normal metal in contact with a superconductor is known to bear superconducting properties 
even if the normal metal has no pairing interaction. The superconducting order in the normal layer 
induced by the proximity effect produces a finite diamagnetic current to expel the magnetic field. 
The quantity of interest is the screening fraction, i.e., how great part of the normal layer contributes 
to the Meissner effect. Theories of diamagnetic response of N-S proximity contact system in both 
the dirty and the clean limit have been already reported.0'i'i'Bi'i'i'& While early experimental 
resultd^lS' were found to be in agreement with the dirty limit theory, some data of more recent 
experiments failed to be described satisfactorily by either the clean or the dirty limit 

theory. 

Belzig et a/J3> proposed a quasi-classical theory of diamagnetic response of a normal layer with 
impurities and pointed out that, even in such a clean system that has a longer mean free path than 
the layer width, temperature dependence of the screening fraction can deviate considerably from 
that of the clean limit theory. Miiller-Allinger et a/.& tried to fit the quasi-classical theory^ with 
their experimental data using the mean free path /n as a fitting parameter. Although they found 
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fairy good agreement in the temperature dependence of the screening fraction with the experimental 
data of the samples with various impurity concentration, there remain some samples which exhibit 
considerably different behavior that cannot be fitted by the theory of ref. |^ Moreover, best fits 
were obtained by using the mean free path a few times smaller than the measured value.!^ 

In the theory by Belzig et al.^ it is assumed that the N-S interface is transparent, i.e., electrons 
pass through the interface without reflection. It is also assumed that the surface scatters electrons 
specularily. The effect of the finite reflection at the interface was examined by Higashitani and 
NagaiS for clean systems and by Hara et al^ for systems with impurities. They showed that the 
finite electron refiection at the N-S interface has strong influence on the screening fraction. 

The purpose of this paper is to study the effect of surface roughness on the screening fraction 
using the quasi-classical Green's function method. The proximity induced order parameter 

extends from the N-S interface into the normal layer with the length scale of the N side coherence 
length (^-^(T). At low temperatures, it will reach the end surface of the normal layer and feel 
the surface roughness. It is expected, therefore, the surface roughness effect is important at low 
temperatures. To focus on the surface roughness effect, we study a clean system and neglect the 
effect of finite refiection at the interface. We show that the surface roughness has a quantitatively 
large effect on the screening fraction. 

This paper is organized as follows. In the next section, we derive a linear response of the 
quasi-classical Green's function to the external vector potential. In §3, we briefiy review the quasi- 
classical theory of rough surface effect by Nagato et al. In §4, explicit form of the quasi-classical 
Green's function in the normal layer of our model system is given. Diamagnetic current and the 
screening fraction are discussed in §5 and 6. The last section is devoted to summary and discussion. 
Throughout this paper, we use the unit h = k-Q = \. 

§2. Quasi- classical Theory of Diamagnetic Current 

In this paper, we consider a clean N-S double layer system as shown in Fig. 1. The end wall at 
z = of the normal layer has roughness, but we assume for simplicity that the N-S interface has a 
translational symmetry in the x-y plane. 

We calculate the linear response of the diamagnetic current to the vector potential using quasi- 
classical Green's function method proposed by Ashida et aZ.lll^ and Nagato et The vector 
potential is taken as A = (yl(z),0, 0), so that the magnetic field is in the y-direction. 



We begin with the Gor'kov Green's function. According to ref. 16, the Gor'kov Green's function 
is given by a product of the slowly varying part and the rapidly oscillating part with the period of 
the Fermi wave length. In the geometry of present interest, the Gor'kov Green's function can be 
expanded in a form 

G(r,r',^„) = E E E G',/5(i^,i^',^,^0e^''-^-''''-^'e'"'^-'''''^', (2.1) 

K,K' a=± /3=± 
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where u;„ is the Matsubara frequency, K is the s = (x, y) component and k = \/2m*E-F — K'^ 
is the 2-component of the Fermi momentum. Note that the slowly varying Green's functions 
Gaf3{K, K' , z, z') are not diagonal in the directional (a) space because of the presence of the surface 
reflection. They are also not diagonal in K space because of the surface roughness. They obey the 
first order differential equationEi'0 

iujn + ^avKP3dz + —A{z)K^ - ) GapiK.K' , z, z) = 5ap6KK'5{z - z), (2.2) 



mc 

where vk = k/m is the z-component of the Fermi velocity, p-^ is a Pauli matrix in the particle- hole 
space and A is the order parameter matrix given by 

A„^f ° ^""^■■^'V (2.3) 

It is convenient to introduce a one-point function Gaf3{K, K' , z) defined by 

Ga/3{K, K', z) ± iaSafsSKK' = -'^^/vKVK'PsGapiK, K' , z±0,z), (2.4) 

where we have taken into account the fact that Ga/3{K, K' , z, z') has a jump at z = z' because of 
the delta function in the right hand side of eq.(2.2). The one-point function obeys an equation 

dzGa/siK, K', z) = i — SKaPzGapiK, K\ z) - i GafsiK, K', z)iK'f3P3, (2.5) 

VK VK' 
6 

EKa = iw„ H A{z)K^ - Aa. (2.6) 

mc 

When K = K' and q = /?, this equation coincides with the Eilenberger equation.© Thus the 
one-point function GapiK, K' , z) can be interpreted as generalized quasi-classical Green's function. 

Since eq.(2.5) is a first-order differential equation, the spatial dependence of GapiK, K' , z) can 
be described by 

G^p{K, K', z) = U^{K, z, z')G^p{K, K' , z')Up{K', z'z), (2.7) 
where Ua{K, z, z') is the evolution operatorEl) that satisfies 

WKdzUa{K,z,z') = -aeKaPzUa{K,z,z') (2.8) 

and the initial condition 

U^{K,z,z) = l. (2.9) 

It is straightforward to calculate from eq.(2.2) the linear response of the Green's function 
6Gai3{K, K' , z, z') to the vector potential A{z): 



6G^p{K,K',z,z')=J2Yl / dz"G^^^{K,K",z,z") 

7 K" "^^ 



-—A{z)K, 
mc 



Gf^{K",K',z\z'),(2.lQ) 



3 



where is the Green's function without magnetic field. In what follows, we omit the superscript 
(0). The diamagnetic current is calculated from the quasi-classical Green's function 6Gaa{K, K' ,z). 
Using eqs.(2.4) and (p?7|) we find& 



6Gaa{K,K',z) 

Uo^iK, z,0)J2Y.\r d^'GW(^, K", 0)4(K", (K" , K' , 0) 

/oo 
dz'Gy [K, K\ Q)L,{K'\ {K", K', 0) 



7 K" 
+ 



Ua{K',0,z), (2.11) 



where 

G'g)(K,i^',0) = Ga(i{K,K',0)±ia6a(i6KK', (2.12) 
CJK,z') = UJK,0, z')—A{z')—p3U^{K, z' ,0). (2.13) 

The quasi-classical Green's function 5Gaa{K, K' , z) is completely determined if we have a knowledge 
of the surface value of Gap{K,K' that can be determined by solving the boundary problem. 

§3. Rough Surface Boundary Condition 

To treat the rough surface effect, we use the random S-matrix theory developed by Nagato 



et alx-M^ Let us consider the scattering of an electron at the Fermi level by the surface. The 
scattering from an initial state {K,—k) to a final state {Q,q) is described by an S'-matrix Skq- 
Since S" is a unitary matrix, we can rewrite S using an Hermite matrix r/ 

S = -\=^. (3.1) 

1 + 17] 

Nagato et a/.llBIll' showed that one can have a formal solution for Gaf3{K, K' ,0) that satisfies the 
boundary condition described by the S'-matrix. The quasi-classical Green's function thus obtained 
is that for a particular configuration of the rough surface, while the surface roughness will be 
distributed randomly enough over the surface. What we are interested in is the quasi-classical 
Green's function averaged over such randomness. 

In what follows, we treat every element of rj as random variable to describe the statistical prop- 
erties of the surface. We assume that rjKQ^s obey the Gaussian statistics with tJkq = and 



(2) r 

VkqVk'Q' = Vkq^Q-k,q'-k'- 

It is instructive to discuss here on the nature of the surface scattering in the normal state. The 
average of the specular reflection amplitude is evaluated within the self-consistent Born approxi- 
mation to 



Skk = - [^r^] =-{l-2r,^ + 2r]^-..-jj,K 

KK 



l + 'vq 

1-ok \ 



(3.2) 
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where obeys an integral equation 



In the same approximation, we can calculate the average of the scattering probability l^j^-Qp: 

\Skq\'^ = \Skk\'^^kq 

(2) 

We adopt a simplest model and put rj^q constant: 



When = 1, it simulates a completely diffusive surface. The specular reflection amplitude Skk 



vanishes and the scattering probability |S'i^Qp is independent of the outgoing momentum Q. Nagato 
et a/.lil'0^ have shown that this model reproduces the diffusive wall boundary condition for the 
Ginzburg-Landau equation of the p-wave state given by Ambegaokar, de Gennes and Rainer.0^ 
The case with W = Q corresponds to the specular surface. Thus, by changing W from = to 
= 1, we can discuss the surface scattering from the specular limit to the diffusive limit. 
Let us consider the random average of 5Gaa{K, K' , z) of eq.(2.11). It is convenient to write 

SGaa{K,K',z) = Uc^iK, z,0)J2Y.\r K", 0)C^{K", z')g[-J {K", K', 0) 

+ dz'Gk~i> {K, K", 0)£^{K", z')G^^J {K", K', 0) 
+ / dz'X^^{K,K",K',z') Ua{K',0,z), (3.7) 

JO 

where Xa^yiK, K", K',z') is defined by 



X^^{K,K",K',z) = (^G^^{K,K",0) - Ga^{K,K",0) 



X C^{K",0,z) i^G^a{K",K',0) - G^a{K",K',0) j .{31 
The average Go,f3{K, K' ,0) in the random S'-matrix model has been given by Nagato et al. 

00) 



G^f,{K,K',0) = Gaf,{K,0)6KK', (3.9) 



G^(3{K, 0) = ( -G^(i^, 0)-i + ^(a - (3) 



+ (g'{K, 0)-i - ia) .^^^^^1 ^ (G'^iK, 0)-i + i/3) , (3.10) 
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where G^{K,0) is the Green's function for the specular surfaceEJ) 
and ax is the surface self energy determined by an integral equation 

The average of Xa-^/iK, K" , K' , z) that corresponds to the vertex correction in the theory of 
impurity scattering can be performedEi^ in the same way as in the calculation of the scattering 
probability of eq.(3.4). 

Y,Y.^c.y{K,K",K',z)=X^{K,z)6KK', (3.12) 

7 K" 

X^iK,z) = (G^(i^,0)'i -ia) |^yA'E^KQnQ(z)y^j (G^(i^,0)-i +ia) , (3.13) 



where 



Yk = , (3.14) 

GS{K,0)-^-aK 

Uk{z)=Yk ^MiK,z)+Y.7jPQUQiz)^ Yk, (3.15) 

MiK, z)=Y^ (g^{K, 0)-i + It) C^{K, 0, z) (G^'iK, 0)-^ - i^) . (3.16) 
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Thus the averaged quasi-classical Green's function becomes diagonal in K space: 



SGaaiK, K', z) = 6Gaa{K, z)6kk'- (3.17) 

§4. Quasi-classical Green's Function in the Normal Layer 

In this section, we discuss the quasi-classical Green's function in the normal layer of N-S proximity 
contact system. In order to focus on the rough surface effect, we assume that the superconducting 
layer is an s-wave semi-infinite superconductor and the normal layer has no pairing interaction, i.e., 

,0 for < z < d, , ^ 

A„(i^,z) = <; " " (4-1) 

L^{z) for a < z < oo. 



We also adopt the simplest model of eq.( |3.6| ) for surface roughness, as a result the surface self 
energy ax is a constant independent of K. 

(2) 

Under these assumptions, one finds from eq.(3.15) that J2q VxQ^Qi'^) ~ because Cj{K,0, z) 
is an odd function of K. It follows that Xa{K, z) = and the third term in eq.(3.7) can be omitted. 
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We further assume that the N-S interface at z = d is transparent. Then, we have for the averaged 
Green's function in the normal layer 



J Z 

roo 

+ dz'G^-^^C^^{K,z')G^+J 



+) 
7a 



where 



^{K,z') = U^{K,0,z')—A{z')^psU}'{K,z',0), 



mc 



VK 



B{K,z') = U^{K,0,d)U^{K,d,z')—A{z')^psU^{K,z',d)U^{K,d,0). 



mc vx 

Here, is the evolution operator in the N (S) layer. 

The explicit form of the Green's function G++ is given 

A 



(4.2) 

(4.3) 
(4.4) 

(4.5) 



G++(i^,0)=i - 



Tr^ 



1 



where 



A = U'^{K,0,d)Ul{K, d,oo)Ul{K, 00, d)U'^{K,d, 0)j 



1 -ifj 



+ icr 



(4.6) 



(4.7) 



and Tr denotes the trace in particle-hole space. We note that the surface self energy ax = o" is 
-ftT-independent in the present model of eq.( |3.6D . The U^{z,z') carries all the information on the 
superconducting layer. It depends in general on the spatial profile of the order parameter A(z). 
But we assume for simplicity that A{z) is a constant and is equal to the bulk BCS value A(T). By 
this simplification, we can have an analytic solution for U^{z,z') as well as for U^{z,z'): 



Uaiz, z') = exp 
Ua{z, z') = exp 
Substituting them into eq.(4.7), we find 



a 
a 



(^nPsiz - z) 
{lOnPZ + Ap2) {z - z] 



(4.8) 
(4.9) 



20 



(17 + ujn)^'^'^ -iA \ 1 - id 

iA {^-ujn)e~'"^ j 1 + 



(4.10) 



where Q, = + A-^ and k = 2uJn/vY^- The surface self energy is determined from eq.(3.11). This 
has to be done numerically for each Matsubara frequency i^n- 
The other Green's functions G-\ , G |- and G are related to (?++ via eq.(3.10). 
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§5. Diamagnetic Current 

In this section, we calculate the diamagnetic current using the quasi-classical Green's function 
derived in the previous section. The diamagnetic current j{z) which flows in the x-direction in the 
present geometry is given by 



i^n K 



= -e(^7V(0)/2)TV / d^sin^ / -typ^Tr psJ^SGaaiK, z), (5.1) 

where A^(0) is the density of states at the Fermi level, is the polar angle of the Fermi momentum 
and vp^ = vp sin 9 cos (p is the x-component of the Fermi velocity. 

As was shown by Zaikini and Higashitani and Nagai,& the diamagnetic current j{z) is constant 
throughout the normal layer and the current-field relation is completely nonlocal. This is charac- 
teristic to the clean system and happens because the normal state evolution operator 11^ {K, z, z') 
commutes with p3, therefore Tr/33 J2a ^^^^{K, z) is independent of z. 

We show below that rough surface effect is large even if the normal layer width d is larger than 
the coherence length ^ and the penetration depth A. In such a case, we can neglect the third term 
in eq.(4.2), because the vector potential A{z) decays exponentially. Now we have an expression for 
the diamagnetic current in the normal layer: 

j{z) = -—n*^A, (5.2) 
mc 

where A is the vector potential averaged over the normal layer and 

TT-N = mdTTT^{{IC/vpJ) (5.3) 

is the effective superfluid density, i.e., the density of proximity-induced superconducting electrons 
in the normal layer. Here, 

/C = iTrp3EG'i+V3G'y, (5.4) 

f^N = Pf/ '^'^'^ is the electron number density in the normal layer and the double bracket means the 
angle average over the polar angle 

Combining eq.(5.2) with the Maxwell equation = —£^d'^A{z), we find the magnetic field B(z) 
in the normal layer when the applied external magnetic field is H: 



r^'^ dOsin^e (•••)• (5.5) 
Jo 



B{z) = H{1- z/X), (5.6) 
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where 



A-^f^ + ll (5.7) 



and Anl is the London penetration depth of the normal layer defined by 



Anl = J^^. (5.8) 

In Fig. 2, we show the surface roughness dependence of n^/nN for (i/,^=20 and 50 = vp/irTc 
the coherence length). The solid curve represents the temperature dependence of n^/riN in the 
diffusive limit and the dashed one represents the specular limitJl* In both the limits shows 
similar temperature dependence but the magnitude is considerably different. As was shown in ref. 
|5|, is exponentially small at high temperatures and begins to grow below temperature ~ (^/(i)Tc 
at which the order parameter extends up to the normal layer end. The magnitude of is strongly 
reduced, in particular, at low temperatures by the diffusive scattering effect. The largest reduction 
is obtained at T = 0. We see that nJ^/riN at T = reaches to unity in the specular case but is 
reduced to about 0.65 in the diffusive limit. We can also see that the suppression rate of n^/riN does 
not depend upon the normal layer width d. This means that, at low temperatures, the rough surface 
effect has to be taken into account even in such thick normal layers that satisfy the conditions d > ^ 
or d > Anl- 

§6. Screening Length 

In this section, we discuss the screening length p defined by 



P = ^£dz{H-B{z)). (6.1) 



This length is equivalent to total diamagnetic susceptibility of the normal layer. The explicit 
expression for the screening length p in clean N-S systems was first given in ref. H. The same result 
holds for the present N-S system and is given by 

3 1 

^/"^ = 4l + 3{X^^/dr{m/n*^y ^^-^^ 

Figure 3 shows the temperature dependence of the screening fraction p/d when d/S, = 20 and 
d/S, = 50. The solid curves represent the diffusive limit and the dashed curves represent the 
specular limit. The temperature dependence of p/d comes only from n^. When n^/riN <C 1 at 
high temperatures, the screening fraction can be approximated by p/d ~ ^{n^ / n^^) / {Xjqi^ / d)"^ and 
is proportional to the effective superfluid density n^. Accordingly, p/d increases exponentially 
with decreasing temperature. At low temperatures, since nJ^/riN becomes of order unity and 
we are considering thick normal layers such that satisfy Anl/c? ^ 1, we may neglect the term 
'i{^NL / d)"^ {n^ / n^) in the denominator of eq.(6.2), so that p/d exhibits saturation to 3/4. These 
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high and low temperature behaviors of p/d are common to the specular and the diffusive limits. 
The difference between the two results is therefore visible only in the intermediate temperature 
range, though the superfluid density nj^ has significant difference at all temperatures. 

Belzig et a/.B!-* discussed the impurity scattering effect on p/d in a same system with our's 
but with specular surface. Hara et a/.@ considered the effect by finite reflection at the interface 
in addition to the impurity effect. The finite reflection has a similar effect to the rough surface 
effect, namely, that dose not alter much the qualitative temperature dependence of p but reduces 
considerably its magnitude. The impurity effect is more involved. In dirty systems, the screening 
fraction p/d increases rather slowly with decreasing temperature (p oc T^^l'^ in the dirty limit) 
and reaches to a larger value than the clean limit one 3/4 at T = 0. The detailed structure of 
the temperature dependence of p is sensitive to the impurity concentration.^ The difference in p 
between clean and dirty systems originates from the fact that the impurity scattering makes the 
magnetic response be locahi'S 

When compared with the results of Belzig et aZ.i'H' and Hara et al.^^ the present analysis suggests 
that the rough surface scattering effect is quantitatively as large as those by impurity scattering 
and by finite reflection at the interface. This means that when a system is comparatively clean 
we have to take into account the rough surface effect simultaneously with those by the impurity 
scattering and by the finite reflection at the interface. 

§7. Summary and Discussion 

We have studied the rough surface effect on the diamagnetic current in the N-S proximity-contact 
system using the quasi-classical theory. The linear response of the quasi-classical Green's function 
to the vector potential has been calculated correctly taking into account the rough surface effect. 

We studied the temperature dependence of the diamagnetic current in the clean normal layer and 
found that the diamagnetic current is significantly suppressed by the diffusive surface scattering 
even in a thick normal layer of which the width is larger than the coherence length ^ = vp/irTc. We 
have shown that in comparatively clean systems the rough surface scattering effect is as significant 
as those by the impurity scattering and by the finite reflection at the interface. Those effects should 
be considered simultaneously for the quantitative comparison with experiments. 

Before concluding this report it is worth mentioning the paramagnetic re-entrant behavior of a 
N-S proximity contact system reported by Mota and co-workers. In the present calculation, 
no such behavior was found. Bruder and Imryllll' suggested that the glancing current along the 
surface of normal metal cylinder contributes to the paramagnetic susceptibility. But the theory was 
criticized for the predicted paramagnetism being too small to explain the experiment, i) Fauchere 
et al. i) pointed out that when the pairing interaction in the normal metal is repulsive, a zero 
energy bound state is formed at the interface and that zero-energy bound state contributes to the 
paramagnetic current.!!^ The bound state energy, however, is not always equal to zero but depends 
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on the reflection coefficient R at the interface. More recently, Maki and HaaaHH^ proposed that the 
normal metals (noble metals) may become p-wave superconductor and generate a counter-current 
in the normal metal. The problem of the paramagnetic re-entrant behavior is still to be studied 
theoretically as well as experimentally. 
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Fig. 1. A geometry of N-S double layer system with a rough surface. The solid and the dashed curves represent the 
order parameter on the S side and on the N side, respectively. 

Fig. 2. The effective superfluid density in the normal layer as a function of temperature T/Tc for layer widths 
d/f = 20 and d/^ = 50. Solid curves represent the diffusive limit. Dashed curves represent the specular limit. 

Fig. 3. Temperature dependence of screening length p for layer widths d/^ = 20 and d/S, = 50. Solid curves 
represent the diffusive limit. Dashed curves represent the specular limit. 
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